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Global Attractivity in Nonlinear Higher Order 
Difference Equations in Banach Algebras 
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("^ ' Abstract 

CN ! 

Non- autonomous, higher order difference equations with hnear argu- 
j^ ' ments of type 



/ ^ ^i^n — i I Qn I / ^ ^i'^n—i 



i=0 \i=0 



^Q , are well-defined on Banach algebras. Their scalar forms with real vari- 

^^ ' ables and parameters have appeared frequently in the literature. By 

generalizing existing results from real numbers to algebras and using 
a new result on reduction of order, new sufficient conditions are ob- 
2 . tained for the convergence to zero of all solutions of nonlinear difference 

equations with linear arguments. Where reduction of order is possible, 
these conditions extend the ranges of parameters for which the origin 
is a global attractor even when all variables and parameters are real 
^ , numbers. 

l> 

(N 

p 

cr^ ■ 1 Introduction 

^^ I Special cases of the following type of higher order difference equation have 

frequently appeared in the literature in different contexts, both pure and 
apphed: 

X- 

?H ■ ^ / '' \ 

.'^: Xn+l = ^aiXn-i + gn{^biXn-i\ , n = 0, 1, 2, . . . (1) 

j=0 \j=0 / 

We assume here that A; is a fixed positive integer and for each n, the 
function (^„ : X — ?> X is defined on a real or complex Banach algebra X with 
identity. The parameters Oj , 6j are fixed elements in X such that 

Ofc / or 6fc / 0. 
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Upon iteration, Equation ([T]) generates a unique sequence of points {a;„} 
in X (its solution) from any given set of k+1 initial values xq, x_i, . . . , x.^ G 
X. The number A; + 1 is the order of the difference equation. 

Special cases of Equation ([1]), in the set of real numbers, appeared in the 
classical economic models of the business cycle in twentieth century in the 
works of Hicks [8], Puu [T7], Samuelson [18] and others; see [21], Section 5.1 
for some background and references. Other special cases of ([1| occurred later 
in mathematical studies of biological models ranging from whale populations 
to neuron activity; see, e.g., Clark [2], Fisher and Goh [5j, Hamaya |7] and 
Section 2.5 in Kocic and Ladas |11] . 

The dynamics of special cases of ([T]) with X = M have been investigated 
by several authors. Hamaya uses Liapunov and semicycle methods in [7] to 
obtain sufficient conditions for the global attractivity of the origin for the 
following special case of ([1]) 

Xn+i = axn + a tanh I 2;n - ^ hxn-i j 

with 0<a<l,a>0 and 6j > 0. These results can also be obtained 
using only the contraction method in [20] • The results in |20j are also used 
in |21j . Section 4. 3D, to prove the global asymptotic stability of the origin 
for an autonomous special case of ([T]) with Oj , 6j > for all i and Qn = 9 for 
all n, where 5 is a continuous, non-negative function. The study of global 
attractivity and stability of fixed points for other special cases of ([1]) appear 
in [6] and {9j; also see [TT], Section 6.9. 

The second-order case {k = 1) has been studied in greater depth. Kent 
and Sedaghat obtain sufficient conditions in [lOj for the boundedness and 
global asymptotic stability of 

Xn+l = CXn + g{Xn - X„_i) (2) 

Also see [22] . In |3] , El-Morshedy improves the convergence results of [10] 
for ([2|) and also gives necessary and sufficient conditions for the occurrence 
of oscillations. The boundedness of solutions of ^ is studied in [19] and 
periodic and monotone solutions of ([2|) are discussed in [23'. Li and Zhang 
study the bifurcations of solutions of ^ in [13]; their results include the 
Neimark-Sacker bifurcation (discrete analog of Hopf ) . 

A more general form of ([2]), i.e., the following equation 

Xn+l = aXn + bXn-1 + S'n(a;n " CX„_i) (3) 



is studied in |24j where sufficient conditions for the occurrence of periodic 
solutions, hmit cycles and chaotic behavior are obtained using reduction 
of order and factorization of the above difference equation into a pair of 
equations of lower order. See [26j for some background on order reduction 
methods. These methods are used in [3] to determine sufficient conditions for 
occurrence of limit cycles and chaos in certain rational difference equations 
of the following type 

aXn + pXn-l + 7 

that are special cases of ^. 

In this paper, by generalizing recent results on reduction of order, to- 
gether with generalizations of some convergence results from the literature 
we obtain sufficient conditions for the global attractivity of the origin for 
([1]) in the context of Banach algebras. These results also extend previously 
known parameter ranges, even in the case of real numbers, i.e., X = M and 
show that convergence may occur in some cases where the functions gn or 
the unfolding map of ([T|) are not contractions. 

Unless otherwise stated, throughout the rest of this paper X will denote 
a real or complex Banach algebra with identity 1 (since there is very little 
likelihood of confusion, 1 also denotes the identity of the underlying field of 
real or complex numbers). For the basics of Banach algebras, see, e.g., [12] or 
|27j . Each Banach algebra is a Banach space together with a multiplication 
operation xy that is associative, distributes over addition and satisfies the 
norm inequality 

\xy\ < \x\\y\ (4) 

with |1| = 1. The multiplication by real or complex numbers (scalars) that 
is inherited from the vector space structure of X is made consistent with the 
main multiplication by assuming that the following equalities hold for all 
scalars a: 

a{xy) = {ax)y = x{ay). 

Elements of type al where a is a real (or complex) number are the 
constants in X. The set M (C) is a real (complex) commutative Banach 
algebra with identity over itself with respect to the ordinary addition and 
multiplication of complex numbers and the absolute value as norm. The set 
C[0, 1] of all continuous real valued functions on the interval [0,1] forms a 
commutative, real Banach algebra relative to the sup, or max, norm. The 
identity element is the constant function x{r) = 1 for all r € [0,1]. The 
other constants in C[0, 1] are just the constant functions on [0,1]. 



An element a; G X is invertible, or a unit, if there is x~^ € X (the inverse 
of x) such that x~^x = 1. The cohection of ah invertible elements of X forms 
a group Q (the group of units) that contains all nonzero constants. For each 
u € ^ if X € X satisfies the inequality 

1 

\x — u\ < 



u'M 



then it can be shown that x £ G. It follows Q is open relative to the metric 
topology of X and contains an open ball of radius l/|u~^| centered about 
each u (z G- Since the zero element is not invertible, ty ^ X. If X is either M 
or C then G = X\{0}. In the algebra C[0, 1] units are functions that do not 
assume the (scalar) value 0. 

2 General results on convergence 

Consider the non- autonomous difference equation 

■^n+l Jn\-^ni ■^n— 1; • • • ) •^n—k) \'^) 

with a given sequence of functions /„ : X^^ — > X. We say that the origin is 
globally exponentially stable if all solutions {xn} of ([5]) in X satisfy the norm 
inequality 

where c € (0, 1) and fi > are real constants such that c is independent of 
the initial values xq, x_i, . . . , x_fc G X. 

The next result, which is true for all Banach spaces (not just algebras) 
generalizes Theorem 3 in |20j . 



Lemma 1 Let X 6e a Banach space and assume that for some real a G (0, 1) 
the functions fn satisfy the norm inequality 

\fn{Co,Ci,---,^k)\ < amaxjl^ol,--- ,|Cfc|} (6) 

for every n and all {^q, . . . ,£,k) S X'^+-'^. Then every solution {xn} of ^ 
with given initial values xq, X-i, . . . , x^^ € X satisfies 

\xn\ < a" ^ max{|xo|, \x-i\, . . . , |a;_fc|}. 
Therefore, the origin is globally exponentially stable. 



Proof. Let fj, = max{|xo|, |x_i|, . . . , |x_fc|}. If {x„} is the solution of ^ 
with the given initial values then we first claim that |a;„| < afi for all n > 1. 
By® 

\xi\ = \foixo,X-i,...,x_k)\ < amaxjlxol,... ,\^~k\} = afi 
and if for any j > 1 it is true that jx„| < afj, for n = 1,2, . . . ,j then 

\xj+i\ = \fj{xj,Xj-.i,...,Xj_k)\ < amax{|xj|,|2;j_i|,..., \xj^k\} 
< a max{/i, O/u} = afj,. 

Therefore, our claim is true by induction. In particular, since < a < 1 
we have shown that \xn\ < a"''^^'^^'fi for n = 1, 2, . . . , A; + 1. Now suppose 
that \xn\ < a^'^^^'^' ^ is true for n < m where m > A; + 1. Then 

|3^m+l I — \Jm\Xmj Xm—lj • • • i ^Jj^—fc j | \ (X maX-[[Xn-t[, |2;jyj_x [,..., \Xm—k | j 

< a/i maxJQ^/^'^+i) , a(™-i)/(^-+i) , . . . , aC'n-fc)/(fc+i) } 

and the proof is complete by induction. ■ 

The above induction argument is used by Berezansky, Braverman and Liz 
in [T] in the case X = R and by Xiao and Yang in [25] in the autonomous 
case {fn = / is independent of n) for general Banach spaces. As we see 
above, this induction argument generalizes to non-autonomous equations in 
Banach spaces. Other approaches that yield convergence results similar to 
Lemma [l] for X = M are discussed by Liz in |14j . 

For X = M Lemma [1] is also implied by Theorem 2 in |15j where Memar- 
bashi uses a contraction argument adapted from Theorem 3 in [_2QJ (ex- 
ponential stability, autonomous case in M). Contraction arguments with 
their geometric flavor are intuitively appealing and they also work for non- 
exponential asymptotic stability; see [20] for the autonomous case and [16j 
which extends the result in [20] to certain non-autonomous equations. 

For a general Banach space the type of convergence is dictated by the 
given norm. For instance, in C[0, 1] with the sup, or max, norm convergence 
to the zero function in Lemma [Tj is uniform. 

Next, define the following sequence of functions on a Banach algebra X 



k / '^ \ 

i=0 \i=0 / 



/n(Co, 6, • • • , Cfc) = >^ aiS,i +gn\}^ biii (7) 



The following corollary of Lemma [T] generalizes previous convergence theo- 
rems proved for the autonomous case with X = M; e.g., the results in [7] or 
Theorem 4.3.9(b) in [2l]. 

Lemma 2 Let (7„ : X — )• X 6e a sequence of functions on a real or complex 
Banach algebra X. Assume that there is a real number a > such that 

\9nm<<y\(\, eex (8) 

for all n and further, for coefficients ai , bi (real or complex) we assume that 

the inequality 

k 

J2(\a^\ + (tM < 1 (9) 

1=0 
holds. Then every solution {x„} of ^^ with initial values xq, X-i, . . . , x_fc G 
X satisfies 

k 

\xn\ < a"'/*^''^"^^ max{|xo|, |2;_i|, . . . , |x_a..|}, a = ^(|aj| + a\bi\). 

i=0 

Proof. If (^0,6, • • • ,6) e X'"'^^ then by the triangle inequality, (j!]) and ([8]) 

k 



k / '' \ 

i=0 \i=0 / 



j=0 

k 



< 



^(|ai| + a\bi 



.i=0 



max{|^o|,---,l6|} 



Therefore, given ([9]), by Lemma [T] the origin is globally asymptotically 
stable. ■ 

Condition ([8]) implies that the origin is a fixed point of ([T]) since it implies 
that 5n(0) = for all n. Except for this restriction, the functions Qn are 
completely arbitrary. Examples of familiar real functions gn that satisfy ([8]) 
include sint, tan~^ t, and tanht; for instance, tanht is used in [7]. Of course, 
gn need not be bounded; e.g., consider t^/(l +t^). 

3 Reduction of order 

Under certain conditions a special, order-reducing change of variables splits 
or factors Equation ([TJ into a triangular system of two equations of lower 
order; see [2^, Theorem 5.6. The next lemma extends that result from fields 
to algebras. 



Lemma 3 Let (?„ : X ^ X 6e a sequence of functions on an algebra X 
with identity (not necessarily normed) over a field T . If for aj,5j € X the 
polynomials 



k 

k—i 



i=0 i=0 

have a common root p ^ Q, the group of units of\, then each solution {xn} 
of (OP in X satisfies 

Xn+l = pXn + tn+1 (10) 

where the sequence {tn} is the unique solution of the equation: 

k-l /k-l \ 

tn+1 = - ^Pitn-i + 5n E 9«*"-« (^^) 

i=0 \i=0 ) 

in X with initial values t^i = x_j — px_j_i G X for i = 0,1, . . . ,k — 1 and 
coefficients 

Pi = p*+^ - aop* tti and qi = 6oP* + &i/0*~^ H Vh 

in X. Conversely, if {t„} is a solution of \11\) with initial values t^i € X 
then the sequence {x„} that it generates in X via U0\) is a solution of ([ip. 

Proof. Define the functions /„ as in ([7]) and for every ^q,vi, . . . ,Vk in X, 
define Co = Co and for j = 1, . . . , A; and fixed j € Q define 

i=l 

Now the change of variables 

tn= Xn- ^Xn-1 (12) 

in Equation ([T|) reduces its order by one if and only if the quantity 

/n(Co,Cl,--- Xk) -7^0 

is independent of ^o ( [25] or [26j, Theorem 5.1). In this case, the above 
quantity defines a sequence of functions 0n(fi, • • • , Vfc) of /c variables that 



yields a difference equation of order k ([25] or [26], Section 5.5). Now, by 
straightforward calculation 



mo, Ci, • • • , a) - 7^0 = (ao - 7)^0 + E "^- (^"')'^o " E^^"' 



-ly--i+l 



^i + 



J=l 



j=l 






(7~'yeo-E(^~')'"^^' 



fi 



ao-7 + X]"j(7 ^) 

J = l 









+1, 



The last expression above is independent of ^o (foi' arbitrary ^o) if and 
only if 7 can be chosen such that 

fc k 

ao-l + Y,aj{j~y = and bo + Y,^j(.^~y = ^■ 
i=i j=i 

Multiplying the two equalities above on the right by j^ yields 

k 
= ao7' - 7'+' + E ^ji^^'y^' = ^(7) 

fc 

so that 7 must be a common root of the polynomials P and Q. 
Now, let 7 = p be a common root of P and Q in Q and define the 
aforementioned functions (/>„ as 

(j)nivi,...,Vk) = fniCoXl,--- ,Ck) - Pia 

hi /hi 

j=l i=l y j=l i=l 

h h / h h 

i=l jf=i \ i=l j=i 



For each i = 1, • • • ,k, since p is a root of the polynomial P it follows 
that 

k 

j=i 

— [P - aoP -•••-fflj-ip )[p ) 

= p' - aQp"'^ ai_i 

Similarly, since p is also a root of the polynomial Q it follows that 

k 

Y, bj{p-'y-'^' = [hp^-^ + fe.+ip^-'-^ + ■ ■ ■ + h,.^p + bk) ip-y-''^' 



[-bop'-hp"-' 6,„i/-+i) (/,-!)* 



= ~bop'-' ~ bip'-^ 6,_i. 

Now, if the quantities pi and qi are defined as in the statement of this 
Lemma then the preceding calculations show that 

k k 

YaJ{P~y-'^'=P^-l and J^ 6,(p-i)^-+i = -g,_i. 

j=i j=i 

Using these quantities the functions </>„ are determined as follows 



(f)n{vi,...,Vk) 



k / '' \ 

i=l \i=l J 



Identifying Vi with the new variables tn-i+i yields a difference equation 
of order k as follows: 

tn+l = (t>n\tni ■ ■ ■ , ^n-fc+l) 



k / k \ 

fe-l /k-1 \ 



i=0 \i=0 / 

which is Equation (jlip . From (J12p we obtain (jlOp and the proof is complete. 
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Remarks. 1. The preceding result shows that Equation ([T]) sphts into 
the equivalent pair of equations (jlOp and (jlip via the change of variables 
(J12p provided that the polynomials P and Q have a common nonzero root 
p in the group of units of X. Equation ([TT]) . which is of the same type as 
(HJ) but with order reduced by one is the factor equation of ([T]). Equation 
(jlOP which bridges the order (or dimension) gap between ([1]) and pip is the 
cofactor equation. 

2. In the special case where 6j = for all i ([1]) reduces to the linear 
non-homogeneous difference equation 



Xn+l 



'^aiXn-i+gn{0) (13) 



with constant coefficients. Since in this case Q is just the zero polynomial, 
Lemma [3] is applicable with p being any root of P in Q. So, as might be 
expected, when X = M the reduction of order is possible if the homogeneous 
part of (|13p has nonzero eigenvalues. The reduction of order of general 
linear equations (non-homogeneous, non-autonomous) on arbitrary fields is 
discussed in [26], Chapter 7. 

3. Solution of polynomials by factorization is problematic in a general 
Banach algebra X due to limitations of the cancellation law. Requiring all 
nonzero elements of X to be units reduces X to either R or C in commutative 
cases and to quaternions in non-commutative cases; see (|27j). For general 
Banach algebras it is possible to determine, for special cases of ([T|), whether 
P and Q have a common root in G', see the two corollaries in the next section. 

4. Equation pT]) preserves another aspect of ([1]). If aj, 6^ G {0} U Q and 
P and Q have a common root p in Q then the numbers Pi,qi in Lemma [3] 
are also units (or zero). 

4 Extending the ranges of parameters 

The solution of Equation (llOp in terms of t„ is 

n 
Xn = p"xo + ^/5"-^t,- (14) 

i=i 

This formula may be used to translate various properties of a solution 
{tn} of ([TT]) into corresponding properties of the solution {x„} of ([T]). This 
is done for equation ([3]) in [23] for X = M. Doing the same for ([TJ more 
generally yields the following natural consequence of combining Lemmas [2] 
andO 
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Theorem 4 Let (^„ : X — t- X 6e a sequence of functions that satisfy (0) for 
each n. Then every solution {xn} of (OP converges to zero if either (a) or 
(h) below is true: 

(a) Inequality (0j holds; 

(b) The polynomials P, Q in Lemma have a common root p & Q such 

that \p\ < 1 and 

fc-i 

Y,i\P^\ + <^\1^\)<^ (15) 

i=0 

with the coefficients pi, qi in the factor equation ill\) . i = 0,1, . . . ,k — 1. 
Proof, (a) Convergence in this case is an immediate consequence of Lemma 

m 

(b) By an appHcation of Lemma [3] we obtain ()lip . Then, given (jlSh . an 
apphcation of Lemma [2] to (|lip imphes that 

\tn\ < a"/('=+iV 

where fj, = max{|io|, |i-i|, • • • , [t-fc+ill with t-i = X-i — px-i-i for i = 

0, 1, . . . , /c — 1 and 

fc-i 

i=0 

with Pi,qi as given in Lemma [3l Since Ip^l < |pp for each j, taking norms 
in ([H]) yields 

n n— 1 / \/(k+l) \ -^ 

\xn\ < iprkoi+^ipr^>,i < \pnxo\+p\prY.\^T^] ■ ^^^^ 

j=l i=o \ 1^1 / 

If ai/{fc+i) ^ IpI then 

lx„!<Ip| |xo| + M«/^ 1^1 K/(^+i)/|p|]-l 

|/9| |Xo|-h/xa a^/{k+i) _ IpI 

Since a, \p\ < 1 it follows that {x„} converges to zero. If a^''^^^^' = \p\ 
then (flBI) reduces to 

I^Jnl < l/^ria^ol + Ail/or^- 

and by L'Hospital's rule {x„} again converges to zero. ■ 

11 



Corollary 5 Let gn he functions on X satisfying ^3^ for all n > 0. Every 
solution of the difference equation 

Xn+l = aXn + gn (boXn + hXn-l H 1" bkXn-k) , (17) 

a eg, hex, bk + 

converges to zero if \a\ < 1 and the following conditions hold: 

boa^ + bia^-^ + b2a''-^ + • • • + 6^ = 0, (18) 

fc-i 
y\boa' + bia'-^ + --- + h\<-. (19) 

i=0 

Proof. For equation (|17p the polynomials P, Q are 

P(^) = e'^+i - ae^ Q(o = &oe' + bi^-' + --- + bk. 



Thus p = a is their common root in Q if (jlSp holds. The numbers Pi,^j 
that define the factor equation (jlip in this case are 

Pi = p*+-^ - ap' = 0, gj = 6oa* + bia'^^ H h 6j 

Thus, if |a| < 1 then by Theorem [5] every solution of (|17p converges to zero. 

■ 

Remark. The parameter range determined by (jlSp is generally distinct 
from that given by ([9]); hence, Theorem H] or Corollary [5] may imply conver- 
gence to when Lemma [2] does not apply and the unfolding map is not a 
contraction. To illustrate, consider the following equation on the set of real 
numbers: 

x„+i = aXn + On tanh {Xn - bXn-k) (20) 

which is a non-autonomous version of a type of equation discussed in [7]. 
Suppose that the sequence {an} of real numbers is bounded by a > and 
it is otherwise arbitrary. Then 

|a„tanht| = |a„||tanht| < lanll^l ^ <^\'t\ 

for all n and ([8]) holds. If < la| < 1 and b = a^ then by Corollary [5] every 
solution of ([20l) converges to zero if 

1 ^, „• l-\a\^ 



j=0 



a ^-^ 
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i.e., if 

On the other hand, applying Lemma [2] to ()20p with b = a^ produces the 
range 

\a\ + a(l + a ) < 1 => ct < p-^ 

1 + \ar 



which is clearly more restricted than the one given by ()2ip . Note that a and 
a may satisfy ([2T]) but with 

|a| +a(l+ |a|^) > 1. 

Corollary 6 Let gn be functions on X satisfying ^ for all n > 0. For the 
difference equation 

Xn+l = aoXn + aiXn-1 H h OkXn-k + gn (Xn " bXn-l) , (22) 

Oi £X,b £ Q, Ofc / 

assume that \b\ < 1 and the following conditions hold: 

ao6^' + ai6^-i + --- + afc = 6^+i, (23) 

fc-i 
^|6^+^-ao6* ai|<l-a. (24) 

1=0 
Then every solution of /i22\) converges to zero. 

Proof. The polynomials P, Q in this case are 

P{0 = C'+' - «oe' ak, QiO =e- b^-'- 

Clearly, Q{b) = and if Equality ([23]) holds then P{b) = too, so 
Theorem m applies. We calculate the coefficients of the factor equation ([TT]) 
as go = 1) % = if i / and 

Pi = b'+^ - aob' Oi. 



Now, inequality (llSp yields (I24p via a straightforward calculation: 

fe-i 
1 > ^ |6*+^ - aob' ai\+a 



i=0 

fe-i 



^|6^+^-ao6^ a^\+a 



i=Q 
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Thus, if |6l < 1 then by Theorem U] every solution of (j22p converges to 
zero. ■ 

As an apphcation of the preceding corollary, consider the case k = 1, 
i.e., the second-order equation 

Xn+l = aoXn + aiXn-l + gn{Xn " bXn-l) (25) 

which is essentially Equation ([3]) on a Banach algebra X. By Corollary [6l 
every solution of (j25]) converges to zero if the functions Qn satisfy dH) and 



beg, \b\<l, aob + ai = b'^, |ao-6I + cr<l. (26) 

On the other hand, according to Lemma [21 every solution of (I25p con- 
verges to zero if the functions Qn satisfy dH) and 

|ao| + Iail+cr(l + [6[) < 1. (27) 



Parameter values that do not satisfy ()27p may satisfy (\26h . For compar- 
ison, if ai = 6^ — aob then (|27|) may be solved for a to obtain 



1 - |ao| - |6||ao - 6| 

'< tth ■ 

This is a stronger constraint on a than a < 1 — |ao — fe| from (|26|). 
especially if b is not near 0. 

Example. Consider the following difference equation on the real Banach 
algebra C[0, 1] of all continuous functions on the interval [0,1] with the sup, 
or max norm: 



ar_ , /3(/3-ar) _^ [\ ( _ J_ 

r + l^"+ (r + l)2 """-^ + 70 '^"V''" ^ + 1' 



Xn+i = ;:7-7-;-X„ + ^^ ^ Xn~l + I 4>n\Xn- -^^-^-^Xn-l ] dv (28) 



where the functions (/)„ : M — > M are integrable and for each n they satisfy 
the absolute value inequality 

|0n(r)| <a|r|, r e^ 

for some a > 0. Assume that the following inequalities hold 

2 + /3-Q 
0</3<l, 3/3<Q<2 + /3, a< ^ 

and define the coefficient functions 

ar /3{l3-ar) /3 

ao('^) = ——r, ai{r) = — — —pr^, b{r' - 



r + 1' ' ' {r + iy ' ' ' r+1 
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Then b G Q, b'^ — aob = l3{/3 — ar)/{r + 1)^ = a\ and the following are 
true about the norms: 

II "^ " 11,1 /^ o 1 

|ao| = sup — — = -, |6| = sup — — - = /3 < 1 

o<r<i r + 1 1 o<r<i r + 1 



and 

ar — /3 



|ao — ^1 = sup 

0<r<l 



r + 1 



max < , p > = < 1 — cr. 



It follows that the conditions in (|26|) are met. Next, since the functions 
Qn : C[0, 1] -^ C[0, 1] in (p5]) are defined as gn{x){r) = Jq (</>„ o x){r)dr for 
r G [0, 1] and all n, their norms satisfy 



/•r pr 

\gn{x)\ < sup / |</>n(x(r))|(ir < sup / a\x{r)\dr < a\x\ sup r 

0<r<lJo 0<r<lJO 0<r<l 



O" X . 



Therefore, Corollary [6] may be applied to conclude that for every pair of 
initial functions Xo,x_i € C[0, 1] the sequence of functions Xn = Xn{r) that 
satisfy ()28p in C[0, 1] converges uniformly to the zero function. In addition, 
it is worth observing that a^ ^ Q and |ao| > 1 if a > 2, in which case (I27p 
does not hold. 

Remark. The preceding corollaries and Theorem S] are broad applica- 
tions of the reduction of order method to very general equations that improve 
the range of parameters compared to Lemma [2j They show that different 
patterns of delays may be translated into algebraic problems about the poly- 
nomials P and Q and their root structures. In some cases a more efficient 
application of Lemma [3] yields a greater amount of information about the 
behavior of solutions than Theorem H] provides. The next result represents 
a deeper use of order reduction in that sense. 



Theorem 7 In Equation I125\) assume that b G G with \b\ < 1 and oq, ai G X 
such that aob-\-ai = 6^. If xq,x-i are given initial values for I125\) for which 
the solution of the first order equation 

tn+l = (ao - b)tn + Qnitn) (29) 

converges to zero with the initial value to = xq — 6x_i then the corresponding 
solution of i25\} converges to zero. In particular, if the origin is a global 
attractor of the solutions of the first order Equation i29\) then it is also a 
global attractor of the solutions of ^25\). 
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Proof. In this case Q{S,) = ^ — 6 so there is only one root b. Now Lemma [3] 
gives Equation (p9]) if P{b) = 0, i.e., if aob + oi = b^. Finahy, we complete 
the proof by arguing similarly to the proof of Theorem l^b), using pi 



Example, consider the following autonomous equation on the real num- 
bers 

Xn+i = axn + b{b - a)xn-i + (Ttanh (x„ - bxn-i) (30) 



where a>0, 0<b<l and a < b. Equation (I29p in this case is 

tn+i = h{tn) h{C) = {a-b)C + a tanh ^ (31) 

The function h has a fixed point at the origin since /i(0) = 0. Further, 
the origin is the unique fixed point of h if |/i(^)| < |^| for ^ 7^ 0. Since h is 
an odd function, it is enough to consider ^ > 0. In this case, /i(^) < ^ if and 
only if (Ttanh^ < (1 — a + 6)^. Since tanh^ < ^ for ^ > we may conclude 
that 

a <l-a + b. (32) 

Given that a < 6, it is possible to choose l<a<l — a + b and extend 
the range of a beyond what is possible with (j26p or (|27p . which require that 
a < 1. In particular, the function utanh^ is not a contraction near the 
origin in this discussion. 

Routine analysis of the properties of h leads to the following bifurcation 
scenario: 



1. Suppose that 6— l<a<6<l and (j32|) holds. Then b — a < 1 and 
all solutions of (13ip . hence, also all solutions of (I30p converge to zero. 



2. Now we fix b, a and reduce the value of a so that a < 6 — 1 < 0. Then 
the function hoh crosses the diagonal at two points r > and — r and 
a 2-cycle {—t,t} emerges for Equation (f3T|) . Note that (j32]) still holds 
when a is reduced, but the origin is no longer globally attracting. The 
cycle {— T, r} is repelling and generates a repelling 2-cycle for (pOj) : see 
|24j or [26], Section 5.5. The emergence of this cycle implies that {t„} 
is unbounded if |to| > t and it converges to if |to| < t. Therefore, 
the corresponding solution {xn} of (pOl) also converges to if 

|xo -bx-il = \to\ < r; 

i.e., if the initial point (x_i,xo) is between the two parallel lines y = 
b^ + T and y = b^ — t in the (^, y) plane. 
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3. Suppose that a continues to decrease. Then the value of r also de- 
creases and reaches zero when 

a = b — a — 1 

i.e., when the slope of h at the origin is —1. Now, the cycle {—t,t} 
collapses into the origin and turns it into a repelling fixed point. In 
this case, all nonzero solutions of (IHTI) and (1301) are unbounded. 



References 

[1] L. Berezansky, E. Braverman, and E. Liz, Sufficient conditions for the 
global exponential stability of nonautonomous higher order difference 
equations, J. Difference Eq. Appl. 11 (2005), pp. 785-798. 

[2] C.W. Clark, A delayed recruitment model of population dynamics with 
an application to baleen whale populations, J. Math. Biol. 3 (1976), 
pp.381-391. 

[3] M. Dehghan, CM. Kent, R. Mazrooei-Sebdani, N. Ortiz and 
H. Sedaghat, Dynamics of rational difference equations containing 
quadratic terms, J. Difference Eq. Appl. 14 (2008), pp. 191-208. 

[4] H.A. El-Morshedy, On the global attractivity and oscillations in a class 
of second- order difference equations from macroeconomics, J. Difference 
Eq. Appl. 17 (2011), pp.1643-1650. 

[5] M.E. Fisher and B.S. Goh, Stability results for delayed recruitment in 
population dynamics, J. Math. Biol. 19 (1984), pp. 147-156. 

[6] I. Gyori, G. Ladas and P.N. Vlahos, Global attractivity in a delayed 
difference equation. Nonlinear Analy. TMA, 17 (1991), pp. 473-479. 

[7] Y. Hamaya, On the asymptotic behavior of solutions of neuronic differ- 
ence equations, Proc. Int'l. Conf. Difference Eq., Special Func. Appl., 
World Scientific, Singapore, 2007, pp. 258-265. 

[8] J.R. Hicks, A Contribution to the Theory of the Trade Cycle, 2nd ed.. 
Clarendon Press, Oxford, 1965. 

[9] G. Karakostas, C.G. Philos and Y.G. Sficas, The dynamics of some 
discrete population models. Nonlinear Anal., 17 (1991), pp. 1069-1084. 

17 



[10] CM. Kent and H. Sedaghat, Global stability and boundedness in Xn^i = 
cxn + f{xn - x„_i), J. Difference Eq. Appl. 10 (2004), pp.1215-1227. 

[11] V. Kocic and G. Ladas, Global Behavior of Nonlinear Difference Equa- 
tions of Higher Order with Applications, Kluwer Academic, Dordrecht, 
1993. 

[12] E. Kreyszig, Introductory Functional Analysis with Applications, Wiley, 
New York, 1978. 

[13] S. Li and W. Zhang, Bifurcations in a second-order difference equation 
from macroeconomics, J. Difference Eq. Appl. 14 (2008), pp. 91-104. 

[14] E. Liz, Stability of non- autonomous difference equations: simple ideas 
leading to useful results, J. Difference Eq. Appl. 17 (2011), pp. 221-234. 

[15] R. Memarbashi, Sufficient conditions for the exponential stability of 
nonautonomous difference equations, Appl. Math. Lett. 21 (2008), 
pp. 232-235. 

[16] R. Memarbashi, On the stability of non- autonomous difference equa- 
tions, J. Difference Eq. Appl. 14 (2008), pp. 301-307. 

[17] T. Puu, Nonlinear Economic Dynamics, 3rd. ed.. Springer, New York, 
1993. 

[18] P.A. Samuelson, Interaction between the multiplier analysis and the 
principle of acceleration. Rev. Econ. Stat. 21 (1939), pp. 75- 78. 

[19] H. Sedaghat, A class of nonlinear second-order difference equations 
from macroeconomics. Nonlinear Analy. TMA, 29 (1997), pp. 593-603. 

[20] H. Sedaghat, Geometric stability conditions for higher order difference 
equations, J. Math. Anal. Appl. 224 (1998), pp.255-272. 

[21] H. Sedaghat, Nonlinear Difference Equations: Theory with Applications 
to Social Science Models, Kluwer Academic, Dordrecht, 2003. 

[22] H. Sedaghat, Global stability of equilibrium in a nonlinear second-order 
difference equation, Int'l. J. Pure Appl. Math. 8 (2003), pp. 209-223. 

[23] H. Sedaghat, On the Equation Xn+i = cXn + f{xn — Xn-i)^ Fields Inst. 
Comm., 42 (2004), pp.323-326. 



18 



[24] H. Sedaghat, Periodic and chaotic behavior in a class of second order 
difference equations, Adv. Stud. Pure Math., 53 (2009), pp. 311-318. 

[25] H. Sedaghat, Semiconjugate factorization and reduction of order in dif- 
ference equations, larXiv:0907r395lV l . 2009 

[26] H. Sedaghat, Form Symmetries and Reduction of Order in Difference 
Equations, Chapman & Hah/CRC Press, Boca Raton, 2011. 

[27] A. Wilansky, Functional Analysis, Blaisdell, New York, 1964. 

[28] H. Xiao and X-S. Yang, Existence and stability of equilibrium points 
in higher order discrete time systems, Far East J. Dyn. Sys. 5 (2003), 
pp.141-147. 



19 



